JOURNAL OF RESEARCH of the National Bureau of Standards - A. Physics and Chemistry 
Vol. 79A, No. 1, January-February 1975 

Equation of State for Thermodynamic Properties of Fluids* 

Robert D. Goodwin** 

Institute for Basic Standards, National Bureau of Standards, Boulder, Colo. 80302 

September 6, 1974 

This equation of state was developed from PVT compressibility data on methane and ethane. The 
highly -constrained form originates on a given liquid-vapor coexistence boundary (described by equations 
lor the vapor pressures and the orthobaric densities). It then requires only five least-squares coeffi- 
cients, and ensures a qualitatively correct behavior of the P (p, 7") surface and of its derivatives, espe- 
cially about the critical point. This nonanalytie equation yields a maximum in the specific heats C,{p, T) 
at the critical point. 
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Symbols and Units 

Subscripts c and _t refer to critical and liquid triple 

points 
Subscript a refers to liquid-vapor coexistence 
a, b, o nonlinear constants in the equation^of state 

B{[)), density-dependent coefficients in the equa- 

C(p) tion of state 

Cor ( T) , molal heat capacity for saturated liquid 
C° V (T) , molal heat capacity for ideal gas states 
C r (p,T) , molal heat capacity at constant density 
C i>(p,T) , molal heat capacity at constant pressure 
d, density, mol/1 

J, the Joule, 1 N-m 

/, the liter, 10 :i m : * 

mol, 16.043 g of CH 4 ; 30.070 g of C 2 H„, (C 12 

scale) 
P, pressure in bars, 1 bar = 10 5 N/m 2 , (1 

atm= 1.01325 bar) 
Pcr(p), the vapor pressure, bar 
R, the gas constant, 8.31434 (J/mol)/K 

p, d\du density reduced at the liquid triple 

point 
o% d/d c , density reduced at the critical point 

r, temperautre, K, (IPTS-1968) 

To~ (p) , temperature at liquid-vapor coexistence 
0(p). defined locus of temperatures, figure 3 

u(p,T), F/7V(p), temperature reduced at 

coexistence 
v, 1/ri, molal volume, 1/mol 

W(p,T), the speed of sound 
x (T) , 777V, temperature reduced at the critical 

point 
x< r (p) , 7V-(p)/7V, reduced temperature at 

coexistence 
Z(P,p, 7 1 ) Pv/RT, the "compressibility factor" 
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1 . Introduction 

A problem of importance for the natural gas in- 
dustry is the prediction of thermodynamic properties 
of Liquefied, multicomponenl mixtures. For the wide 
range of compositions encountered, it may be neces- 
sary to utilize accurate properties of the pure com- 
ponents. We therefore believe that it will be helpful 
to have a relatively simple and rational equation of 
state of identical form for each component, such as 
the equation of state described below. 

This equation originates on a given, liquid-vapor 
coexistence boundary, thus eliminating the long- 
standing problem of consistency between equations 
of state and this independently-derived envelope 
[26]. It ensures a maximum in the specific heats 
C,(p,T) approaching the critical point, qualita- 
tively consistent with experimental behavior near, but 
not necessarily exactly at this pole [16|, and it has 
only five arbitrary, least-squares coefficients. Experi- 
mental compressibility data for methane and ethane 
have been used at densities to the triple-point liquid 
density, temperatures to twice the critical, and 
pressures to 350 bar or greater. 

Our objective in the present report is to give a con- 
cise description of this new type of equation of state 
which, with full documentation, has been presented 
in previous publications [14, 15]. 

For methane we have shown in [14| that an equation 
similar to that described below yields calculated 
specific heats, C (r (T), C v {p,T), C p {p,T), and speeds 
of sound in acceptable agreement with experimental 
data, without any weighting of the equation of state 
to those data. In a current report on ethane [15] 
we compute provisional thermodynamic and related 
properties by means of the simpler equation of state 
(1) described herein. 

Symbols and units used here are given in a List. 
Fixed-point values from [14, 15] are given in table 1. 
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For equation of state (1) the gas constant is This monotonic behavior (nonnegative slopes) is 
R =(0.0831434)* (d t ) bar/K, consistent with use of difficult to achieve from equations of state, yet is very 
the dimensionless density, p = d\d t . important for such computations as — 

C P (P, T)=C V (P, T) + T-(dPldT)*l(dPldp)lp*, (b) 

Table 1. Fixed-point values from [14, 15] 

W( P ,T) = [C„-(dPldp)IC r ]W (c) 



Triple point 



Density mofy 



Vapor 

Liquid 

Temperature K. 

Pressure bar a . 



Critical point 

Density mol/1. 

Temperature K. 

Pressure bar. 



Metha 



1.567 865-10- 
28.147 
90.680 

0.1174 35675 



10.0 
190.555 
45.988 



Ethane 



1.35114-10 6 
21.68 
89.899 

1.009 906-10" 5 



6.74 
305.37 
48.755 



a Precision required for vapor-pressure equation. 

Various methods for utilizing the equations of state 
of pure components to derive properties of mixtures 
are described in a number of recent publications, e.g., 
[1, 4a, 17, 18, 19, 20, 21, 21a, 22, 23, 24, 271. 




T c T — 

FIGURE 1. The locus of isochore inflection points. 

2. The Equation of State 

The P(p, T) surface and equations of state are 
described in several reports, e.g., [6, 7, 10, 12, 14, 15, 
20, 24, 28]. Figure 1 shows the qualitative behavior of 
isochores as indicated by Rowlinson [20], needed to 
give a calculated maximum in C v (p, T) at the critical 
point via the isothermal computation — 



AC r = -T- J (d*PldT 2 )-dplp 2 . 



(a) 



Figure 2 shows the well-known zero slope and 
curvature of the critical isotherm at the critical point. 




FIGURE 2. Behavior of the critical isotherm. 

The liquid-vapor envelope, T (T (p), figure 3, is an 
important boundary of the P(p, T) surface for the equa- 
tion of state. We constrain the equation to this bound- 
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Figure 3. Behavior of the locus 6(p). 

Point C is the critical point, and T„(p) is the liquid- vapor coexistence envelope. 
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ary as follows, by use of the vapor pressures and the 
coexistence temperatures T (r (p) formulated in the 
appendix. 

For any density (isochore), obtain the coexistence 
temperature from the function T (T (p). Use this to 
obtain the vapor pressure P (r [T a (p)] as a function of 
density, thus defining the equation of state at coexist- 
ence. By subtraction one then obtains an equation of 
the type of eq (1). The melting line is not a part of this 
equation of state 

P-P a (p)=pR'[T-T.(p)]+p*RT,-F(p,T), 

where, 

F( P , T) =B(p)-4>( P , T) + C(p)-V(p, T). (1) 

The functions B (p) and C)p) are polynomial coeffi- 
cients to be found by least squares. The temperature- 
dependent functions ^(p, T) and ^ (p, T) must be 
defined to be zero on the coexistence boundary at 
T=T <T (p) , as shown by figures 4 and 5. 




T — 



Figure 4. Behavior of the function <J>(p, T). 
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Figure 5. Behavior of the function ^(p, T). 
Equation (1) has, in addition to pRT, only two tem- 



describe the sigmoid shape of isochores in the density 
range p c < p < 2-p r , figure 1. 
The first of these functions is — 



4>(p, T)=x l ' 2 -ln[u(p, T)] 9 



(2) 



where a(p, T) =TIT tr (p). It is linear (d 2 &ldT 2 = 0) 
everywhere on the coexistence boundary at u=l, 
figure 4. It therefore gives a critical isochore which is 
linear at the critical point because C(p) = by defini- 
tion along this isochore, eq (6). 

The second of these functions, ^(p,! 7 ), is defined 
such that eq (1) will yield a maximum in the specific 
heats C r (p,T) at the critical point via eq (a). We 
first define the arguments — 

o>(p,T) -6- [7Mp)-l]. 

av(p) -8- [T, r (p)lO(p)-l], 

where 8 is an arbitrary constant, and 0(p) is our 
locus of temperatures inside the coexistence en- 
velope, figure 3, 



6(p) =7V(p)-exp[-a-/(o-)]. 



(3) 



The function f((r) here is normalized to unity at the 
liquid triple-point density — 

f(a) -|<r-l|*/(«r,-l)», 

where cr, = d//d c is a constant. 
Function ^(p, T) now is defined as the difference, 



V(p,T)=*(p t T)-*„(p), 



(4) 



such that ^ = 0at coexistence, T=T (T (p). Component 
functions, i//(p, T), are designed to give infinite curva- 
ture (d 2 \fj/dT 2 ) at the origin, co = 0, 

i//,(p,r) s [l-o) -lnU + l/o))]/*, (4-a) 

i/Up) = [1 - oj <T • ln(l + Ha) (r )]/x <r . (4-b) 

Figure 5 shows behavior of ^(p.T). Sufficiently far 
away from the critical point it behaves roughly like 
l/T 2 , found in the well-known, Beattie-Bridgeman 
equation of state. 

Behavior of the coefficients B(p) and C(p) in eq 
(1) is shown by figure 6 for methane. The following 
polynomial representations have been developed 
tediously by trial, 

B(p) m Bo + fli ' P + B 2 -p 2 /(l + & -p 2 ), (5) 

C(p) - (o— 1) -(cr- Co)-(Ci + C,-p). (6) 



The sign of the curvature of isochores (d 2 P/dT 2 ) 
perature-dependent functions, which is the minimum at the coexistence boundary is determined uniquely 
number of functions (monotonic in T) needed to by the sign of C(p), figure 6, because <J>(p,r) is linear 
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FIGURE 6. Behavior of the coefficients B(p), C(p) for methane. 

on this boundary. The root in C(p) at cr = 1.9 was found 
by least squares for methane and for ethane. It then 
was introduced as the constant Co in eq (6). This con- 
straint is valuable because, under various conditions, 
we often have failed to obtain any such root from 
PVT data by least squares. 
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FIGURE 7. Presumed behavior of C(p) for hydrogen, reflecting 
observed positive isochore curvatures in compressed liquid states at 
the lowest temperatures [7]. 



Figure 7 shows the presumed behavior of C(p) for 
hydrogen, needed to give the observed positive curva- 
ture of isochores in compressed liquid states at the 
lowest temperatures [3, 7]. The behavior of figure 7 
is described by a double root at C , e.g., 

C( P )=C l ■ (cr-1) • ((t-Co) 2 , 

or without a root by — 

C{p) = (cr-1) • (d + C 2 ' p + C 3 -p 2 ). 

Table 2 presents constants of eq (1) for methane and 
ethane. Individual deviations of experimental PVT 
data from eq (1) are given for ethane in [15]. For 
methane with eq (1) they are comparable with those 
obtained in [14| using nine least-squares coefficients. 
Complete computer programs are reported in [15] 
for adjusting the vapor-pressure equation, the ortho- 
baric-densities equations, the coexistence-tempera- 
tures equations, and the equation of state (1). 



Table 2. 


Constants for equat 


ion of state (1 ) 




Methane 


Ethane 


a 
b 
5 

Co 




2 
1 
1/2 
1.9 


2 
1 
1/2 
1.9 


B 

B x 
B 2 




1.5082 12989 
0.6544 90304 
4.1320 82291 


1.8481 67996 
1.5697 04511 
5.5601 86452 


c> 




-0.7654 09076 
-0.0590 88717 


-1.0428 42462 
+ 0.2249 78299 


N 


% 


756 
0.42 


562 
0.57 



3. Comments on the Equation 

Behavior of the critical isotherm from eq (1) at the 
critical point is deduced as follows. The functions 
P (r [Tcr(p)] and <I> (p,T) depend directly upon T (T (p), 
which gives the overpowering factor exp [— y\ \ a — 1 1 ] 
for derivatives with respect to density at the critical 
point, o*=l. The function ^(p, T) has a finite third 
derivative because it depends also on 0(p), eq (3). 
Its coefficient, C(p), however, is zero at the critical 
density, eq (6). The first, second, and third derivatives 
of eq (1) therefore are zero at the critical point. 

Detailed examinations of this isotherm from eq (1) 
show, however, that small variations in the assigned 
(p c , T c ) critical point give small irregularities (nega- 
tive slopes) nearby at a ^ 1. We find that, given an 
accurate value of T c , eq (1) serves to find the critical 
density which yields a well-behaved critical isotherm 
[15]. For methane and ethane the value of the critical 
density obtained by this method is roughly 1 percent 
lower than estimated by the conventional procedure 
of extrapolating the rectilinear diameter to the critical 
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temperature [14, 15, 20], but in each case falls within 
bounds of uncertainty in published works. 

Ioschore inflection points figure 1, calculated from 
eq (1), are obtained as the difference of second deriv- 
atives (versus T) from the functions composing 
F(p, T). We therefore expect high sensitivity to the 
analytical forms of 3>(p, T) and ^(p, T). Variations 
of these forms might improve accuracy in representing 
PVT data. In the following we describe two alternative 
functions for ^(p, T) from among many different 
functions investigated both for <£(p, T) and for 
^(p, T). We then compare computed specific heats 
obtained via these two forms in the modified equation 
of state described below. Each of the following com- 
ponent functions must be differenced, as indicated by 
eq (4), to obtain ^(p, T) for the equation of state. The 
argument for each of the following functions is changed 
to — 

a)(p,r) = [l-0(p)/7]. 

In (7) the adjustable exponent is 1 < e <2, and the 
adjustable coefficient was selected by trial to be 8= 1/e, 



<Mp, T) - 1- (w-o-w<)/(l-o). 



(7) 



This function approaches zero at high temperatures 
in proportion to 1/T 2 , as seen by expanding co e . By 
trial with methane and ethane PVT data, we found 
e = 3/2, and hence d^/dT 2 behaves like (l/w) 1/2 on 
approach to the origin, co = 0. 

In (8) we use the arbitrary constant ^ 8 ^ 1 to give 
relative weighting to two terms behaving, respectively, 
like 1/7 and \/T 2 at high temperatures, 

i// 3 (p, T) =8-(0/r) + (l-8)-[l-G> + o>-ln (o))]. 

(8) 

The last term in (8) was explored for hydrogen [10]. 
Equation (8) can be simplified for computations. This 
function gives d 2 i// 3 /oT 2 ~- (1/w) on approach to the 
origin, cu = 0, similar to definition (4-a) above. 

With above functions, we find that coefficients B\ 
and C2 of eq (1) become nonsignificant for methane, 
leaving an equation of state with only three arbitrary, 
least-squares coefficients. 

For computations on mixtures, the simplest possible 
equation of state is desirable. We have modified (1) 
such that it is suitable for corresponding-states adapta- 
tions, by specifying that densities shall be reduced at 
the critical point (hence R = 0.0831434 -d,., bar/K), 
and by introducing the critical pressure. The functions 
<I>(p, 7") and ^(p, T) remain unchanged — 

P-PAp) = <rR-[T-T (r (p)]+P<-<T 2 -F(p, T), 

(1-A) 



B(o-)=£ ( , + £,-o- 2 /(l + &-o- 2 ), 



C(<r)"C,-(o— l)-(o— Co). 



(5-A) 



(6-A) 



Table 3 for methane gives the rms of relative density 
deviations for authors in [14], corresponding to each 
function i//(p, T) described above and, on the bottom 
line, the mean of combined pressure deviations. 
Tables 4 and 5 give the constants for eq (1-A), so that 
it will be possible to compute specific heats. 

TABLE 3. Methane density deviations, tins percent 



Equation of state 


eq(l) 


eq ( 


-A} 






Function i//(p, T) 


i//,. eq (4-a) 


1//2, eq (7) 


>//:<. eq (8) 




Authors 


N 


0.06 
.06 
.51 

.42 


0.12 
.13 
.60 

.46 




Virial eq, [14] 

Douslin et al. [4| 

NBS [141 

Mean AP/P, %... 


46 
171 
539 

7S6 


0.16 

.10 
.50 

.45 



TABLE 4. Constants for equation of state {I -A), using i//»(p,' 

eq (7) 





Methane 


Ethane 


b 

a 


e 


1/8 

1/2 
2/3 
3/2 


1/8 

1/2 
2/3 

3/2 


Bo 

B 2 


1.9894 21671 
0.7924 35706 


2.2373 56347 
0.9304 97491 


Co 

c 1 


1.9 
-0.8309 40825 


1.9 
-1.0494 11810 


N 


756 
0.46 


562 
0.79 



TABLE 5. Constants for equation of state (/—/), using i//$(p,T) 

eq (8) 





Methane 


Ethane 


b 

a 
b 


1/8 

1/2 
1/3 


1/8 

2/3 
1/2 


Bo 
Bi 


1.9794 55787 
0.8017 24817 


2.1805 54662 
0.9814 53571 


Co 

c x 


1.9 
-0.9465 70360 


1.9 
-1.3543 33999 


N 
AP/P, % 


756 
0.45 


562 
0.75 



The interaction specific heats (C r — CJ?) for methane 
have been compared along the critical isochore, as com- 
puted by equations (a) and (1-A), using the functions 
i// 2 (p, T) and i//.j(p, T). Experimental behavior of spe- 
cific heats suggests that these data should be described 
very near the critical temperature by — 



(C r -C°)~ (TIT c -l)-\ 



(9) 



75 



where exponent < n < l,and usually, 0.05 ^ n ^ 0.15, 
[16]. 

Figure 8 shows the computed results via i//2(p, ^K 
eq (7). Their qualitative behavior is correct, but they 
cannot be represented by (9). Instead, they are 
described accurately to 210 K by — 

(C v -C^ = 15.1- 18.0- {TIT c -l)V\ JlmolIK, 
suggesting that they might become finite at the critical 
temperature. 



Using i//3(p, 7 1 ), eq (8), on the other hand, yields 
results on the logarithmic plot of figure 9. The slope, 
n ~ 1/3, is close to that observed for our experimental 
oxygen specific heats [8]. 

We conclude that i//2(p, T) gives too weak a curva- 
ture to isochores from (1-A) approaching coexistence, 
and that the forms i//i(p, 7) and ifj 3 (p,T) may be pref- 
erable, despite the logarithmic infinity in their first 
derivatives versus T at the origin, co = 0. 
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FIGURE 8. Computed interaction specific heats (C v — C°) of methane 
along the critical isochore, via eq (1-A) and \j/ 2 (p, T). 
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4. Conclusions 

As argued above, the present equation of state 
(1) is rational because we understand the purpose of 
most of its component parts. Inclusion of the vapor 
pressures and orthobaric densities in this equation 
merely incorporates these physical properties which 
almost invariably must be used for a consistent net- 
work of thermodynamic functions, and eliminates 
the long-standing problem of continuity at the co- 
existence boundary [26]. 

We believe eq (1) to be valuable for thermal computa- 
tions because it ensures an inherently correct be- 
havior of the P(p,T) surface, giving a maximum in the 
specific heats C v (p,T) at the critical point. The very 
small number of arbitrary least-squares coefficients 
(five, and possibly only three) facilitates compari- 
sons of the equation of state for different substances, 
and may be attractive for work on mixtures. 

We cannot expect eq (1) to represent some high- 
precision PVT data as well as equations with a much 
greater number of arbitrary constants. As the preci- 
sion of "good" PVT data probably is often much bet- 
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ter than absolute accuracy, especially in the critical 
region, it would appear to be a self-defeating exercise 
to strive for the ultimate "representation" while ignor- 
ing essential features of the P(p,T) surface. Any in- 
accuracies in the given, liquid-vapor, P — p — T 
boundary for eq (1), however, will be propagated along 
calculated isochores because eq (1) originates on this 
boundary. 

Equation (1) almost certainly cannot be integrated 
analytically to express derived properties in closed 
form. It therefore would not be convenient for multi- 
property analysis [2|. Equations amenable to in- 
tegration, however, probably are not accurate in the 
critical region, which influences a large fraction of 
the P(p,T) surface. To some extent, this must create 
a need for multiproperty analysis. 

For methane with an equation of state similar to 
eq (1), on the other hand, we have used a minimum of 
specific heat data only to compute around the critical 
point for high densities near the critical temperature 
(fig. 3). In all other regions we computed specific heats 
and speeds of sound "a priori" via ideal gas specific 
heats and the equation of state, finding acceptable 
agreement with experimental data. We concluded 
that if such data do not exist experimentally, they could 
be estimated via the present type of equation of state 
for many but probably not all regions of the P(p,T) 
surface [14]. 

5. Appendices 
5.1. The Vapor Pressure Equation 

The original form of our vapor pressure function [9] 
is satisfactory for methane [14], but for ethane it has 
been necessary to add the term d-x 4 [15]. Define the 
argument- 
s') = (l-TtlT)l(l-TtlTc), 
when the function is — 

\n{PIP t ) = a-x+b-x 2 + cx 3 + d'X 4 + e-x'(l-x) 3 l 2 . 

(5.1) 
Table 6 presents the coefficients. 

TABLE 6. Constants for vapor-pressure equation (5.1) 



Methane 



4.7774 8580 
1.7606 5363 
-0.5678 8894 
0.0 
1.3278 6231 



Ethane 



10.7954 9166 

8.3589 9001 

-3.1149 0770 

-0.6496 9799 

6.0734 9549 



5.2. The Orthobaric Densities 

The following expressions are constrained to the 
boundaries, the triple- and critical points [15]. In each 
case the basic behavior is described by — 

Y (p, T) = const., 



and polynomials are selected to represent small devia- 
tions. For the saturated liquid, define the variables, 

x(T) = (T c -T)l(T c -T t ), 

y(p) = (p-pc)l(pt-pc), 

Y(p,T)-(y-x)K*-x), 

when the function is — 

Y=a + b-x 2 l* + c-x. (5.2) 

For the saturated vapor, define the variables, 

x(T) - (T c /T-l)l(T c IT t -l), 

yip) = In (p el p)lln(pclpt) 

Y(p,T) - (y-*)/(*«-«), 

when the function is — 



Y=A 1 +'2 A, ■ x'l*. 

i = 2 

Table 7 presents constants for (5.2) and (5.3). 

TABLE 7. Constants for orthobaric densities equations 



(5.3) 



Metha 



Ethi 



Saturated liquid densities, eq (5.2) 



e 


0.36 


0.33 


a 


.8595 3758 


.7219 0944 


b 


.0243 6448 


.2965 7790 


c 


-.0268 5285 


-.3003 6548 





Saturated vapor densities, eq (5.3) 


€ 


0.41 


0.39 


Ay 


.4171 4211 


.2158 7515 


A 2 


-.5194 9762 


-.0852 2342 


A, 


1.2077 7553 


-.6152 3457 


A 4 


-1.4613 0509 


.2545 2490 


A-, 


0.5765 8540 


.1517 7230 



5.3. The Coexistence Envelope, l (r (p) 

This envelope is shown by figure 3. For the equation 
of state (1) we obtain the coexistence temperature for 
any density from the following analytical descrip- 
tions. If the coexistence density must be found at a 
given temperature, it is obtained (for consistency in 
the equation of state) from the following expressions 
by iteration, using eqs (5.2) and (5.3) only to estimate 
an initial density. 

An important feature of the present formulations 
[15] is that derivatives — 

d"Tcj/dp" 

of all orders, n, are zero at cr = l, the critical point. 
We describe T a {p) in two parts, according as a ^ 1. 
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This simplifies constraint to the boundaries (liquid and 
vapor triple points). For each range the dependent 
variable is — 

Y[T, T (p)] '(TJTa-DKTJTt-l), 

and we use the following function, infinite at the critical 
density 



U(<r) 



- r [l/|o— l|-l/|crr-l|] 



where crt=dtldc is a constant, and dt refers to vapor 
or liquid at the triple point according as a ^ 1. 

For the liquid range at cr ^ 1 the equation is — 

\ n (Y) = U(a) + j^A i - (o-'-o-p. (5.4) 

/= 1 

For extremely low densities in the vapor range at 
o"^ 1 we modify the above expression as follows. 
Define the variable 

r(a)=ln(l+s/o-)/ln(l+s/o-r), 

where s is an arbitrary constant. Our equation for the 
vapor range now is — 

ln(F) = I/(cr)+flo-ln(r)+fli- (cr 1 ' 3 -^, 1 ' 3 ) 

+ B 2 • (cr 2 ' 3 - of 3 ) + J Bi • (o-*- 2 - o-;- 2 ) . 
; = 3 

(5.5) 
Table 8 presents constants for (5.4) and (5.5) 



Table 8. Constants for the To-(p) equations 



Methane 



Ethane 



Saturated liquid, eq (5.4) 



y 


1/2 


1/2 


A, 


11.4317 7230 


23.7245 1840 


A, 


-3.8765 9480 


-14.8860 5161 


A, 


0.5378 8326 


5.4317 7443 


A 4 


.0 


-1.0715 0566 


A 5 


.0 


0.0913 5183 



Saturated vapor, eq (5.5) 



s 

B, 

B l 

B 2 

B, 


1/4 
0.9034 9557 
.0 
.0 
-.3834 4338 


1/4 

0.8681 0517 

.0151 6978 

-.7296 0432 

1.0096 5493 


B 4 
B, 
B« 
B 7 


-3.9210 8638 
6.2600 3837 

-9.3296 0083 
5.6060 2816 


-8.7340 2710 

21.1071 2823 

-31.4499 4087 

17.8637 0397 



Anneke L. Sengers emphasized to us some time ago 
the importance of the critical region for the entire 
equation of state, thus motivating present develop- 
ments. In this laboratory, R. D. McCarty provided the 
essential least-squares program, and we are indebted 
to D. E. Diller and L. A. Weber for discussions and 
valuable suggestions. The American Gas Association 
very kindly has supported this work. 
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